Order-disorder transition in the Cd-Ca cubic approximant 
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ABSTRACT 

Recent experiments discovered an order-disorder transition occuring at low 
temperatures in large unit 1/1 cell cubic approximants of the stable Cd-based binary alloy 
quasicrystals. The transition is related to correlations among orientational degrees of 
freedom whose separations are around 12 A. We analyze the interactions between the 
degrees of freedom using ab-initio calculations for Cd-Ca alloys and derive an equivalent 
antiferromagnetic Ising model which shows a similar phase transition. However, the 
calculated transition temperature is higher than observed experimentally, indicating that 
the actual structure and its order-disorder transition are more complex than originally 
proposed. A side-benefit of our study is the discovery of a canonical-cell decoration model 
for the Cd-Ca icosahedral phase. 

INTRODUCTION 

Stable icosahedral quasicrystals occur in the compounds Cds.rYb and Cds.yCa [1, 2, 3]. 
In each case the phase diagram of the binary alloy contains a 1/1 cubic approximant at a 
close-by composition [4, 5]. The CdgYb structure [6] can be represented in a conventional 
simple cubic unit cell with a lattice parameter of a = 15.7 A. The Pearson symbol for 
CdgYb is cI176 indicating that it is cubic, body-centered, with 176 atomic positions per 
simple cubic unit cell. Since 176 is not a multiple of the basic 7-atom stoichiometric unit (6 
Cd and 1 Yb), there must be partial occupancy. In fact, the unit cell contains 144 Cd 
atoms, 24 Yb atoms and 8 vacancies. Sets of 4 vacancies alternate with sets of 4 Cd atoms 
among vertices of a cube to form tetrahedra located at the center of a dodecahedral cluster 
(see Fig. 1). There are two such clusters per simple cubic cell, for a total of 8 vacancies. 






Figure 1: (Left) Innermost Cd2o dodecahedral shell centered by Cd4 tetrahedron, 
viewed along cubic [100] axis. (Right) Two alternate tetrahedral orientations (dark=Cd, 
light = vacancy) according to cI176 structure [6]. 



The proper structure type of CdeCa has recently been debated. Initial reports [7] 
assigned it to prototype CdeY with Pearson symbol cI184. This structure is equivalent to 
the cI176 structure of CdgYb except for the central Cd4 tetrahedra. In cI184 these Cd 
atoms occupy 4 out of 12 vertices of a cuboctahedron (resulting in 6 possible oricntatons) 
instead of 4 out of 8 vertices of a cube (resulting in 2 possible orientations). After the 
discovery of stable Cd-Ca quasicrystals some researchers [8, 9, 10] suggested the structure 
type of CdgCa might actually be CdgYb.cIlTG. The full story has yet to be resolved. 
Indeed a recent study of MCde for a variety of metal atoms M finds evidence for yet 
additional site types [11]. 

In order for the overall structure to be truly body-centered, it is necessary that either 
the tetrahedra at each center be identically aligned or else that all body centers be 
randomly oriented. In each case a body-center translation leaves the structure invariant. 
However, recent experiments indicate order-disorder transitions in which a supercell 
ordering develops at low temperatures, breaking the body-centered cubic 
symmetry [10, 12, 13]. 

Estimated transition entropies close to A;^ In 2 suggest that the cluster center 
tetrahedra have two equivalent orientations. Both should occur randomly at high 
temperatures, and freeze into some definite ordered pattern at low temperatures. Such a 
local two-state system can be modeled using the Ising model, with the two spin states 
representing the two cluster orientations. Because the cI176 structure of CdeYb possesses 
precisely two tetrahedron orientations, we began our initial study of the order-disorder 
transition using this structure. 

TOTAL ENERGY CALCULATIONS 

We carry out ab-initio calculations using the plane-wave program VASP [14, 15] which 
yields reasonably accurate total energies. This approach uses ultrasoft 
pscudopotcntials [16] or PAW potentials [17] to represent the effective interaction of 
valence electrons with ionic cores, and solves the many-body quantum mechanical band 
structure of these electrons using electronic density functional theory. We choose to model 
Cd-Ca rather than Cd-Yb because the alkali earth element Ca is easier to treat from first 
principles than the rare earth element Yb. 

First, we use VASP to reproduce the sequence of low-temperature stable phases in the 
estabhshed Ca-Cd binary phase diagram [5] . To do this we calculate the cohesive energy 
for each known structure, and several hypothetical ones. Each structure is fully relaxed in 
both imit cell parameters and atomic coordinates. All energies are converged to an 
accuracy of 1 meV/atom or better by increasing the k-point mesh density. These 
calculations use PAW potentials in the generalized gradient approximation and a constant 
plane- wave energy cutoff of 274 eV. 

Subtracting each cohesive energy from the tie-line joining the pure elements in their 
ground states yields enthalpies of formation (at T=OK). Enthalpies of all known and many 
hypothetical Ca-Cd structures are plotted in Fig. 2. We label each structure with its name 
followed by its Pearson symbol. 

Agreement between our calculation and the established phase diagram requires that all 
known low temperature structures lie on the convex hull of enthalpy versus composition. 
Additionally, all hypothetical structures must lie above the convex hull, as must all known 
high temperature, high pressure and metastable phases. Agreement is perfect, except for 
one seeming difficulty with CaCd2 which we now address. 

The low and high temperature CaCd2 phases are reversed in energy relative to the 
experimental report [5]. However, the transition between the hP12 and oI12 variants has 
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Figure 2: Enthalpies of formation of Ca-Cd compounds. Notation: heavy circles indicate 
known low temperature phases; light circles indicate known high temperature phases; squares 
indicate either structures not reported in Ca-Cd system. Line sigments connect vertices of 
convex hull. Legends hst calculated stable structures (left) and unstable structures (right). 
Inset shows details near quasicrystal-forming composition. 

not been well established. Most likely, according to our findings, the presumed low 
temperature hP12 phase is actually a metastable phase, and the nominal high temperature 
oI12 phase is actually stable all the way to low temperatures. 

The established phase diagram [5] lists a phase CaaCdiy, of unknown structure, which 
is claimed to exist from T=OK up to melting, at a composition extremely close to CaCdg. 
Thermodynamic rules governing alloy phase diagrams suggest that either CdeCa or this 
unknown phase should actually be stable at high temperature only, because the probability 
is low that distinct phases coexist over an extended temperature range when they are close 
in composition. 

One candidate for the Ca3Cdi7 phase (based on atomic size ratio and chemical 
similarity) is Sc3Zni7.cll60. This structure differs from CaCdg.cIlTG because the innermost 
Cd4 tetrahedra (see Fig. 1) are missing. We find the energy of this structure is substantially 
above the convex hull, confirming that occupancy of the inner shell is energetically 
favorable but leaving open the question if occupancy by four Cd atoms is the optimum. 

Other candidates for the structure of the phase named Ca3Cdi7 are the cluster of 
additional structures just to the left of CdgCa in Fig 2. These are postulated icosahedral 
quasicrystal approximant structures based on canonical cell decorations, as discussed 
below. The prohferation of many distinct nearly degenerate structures, at energies slightly 
above the convex hull, is consistent with the proposal of entropic stabilization of 
quasicrystals. We conclude that the observed phase named CasCdiy is most likely the 
icosahedral phase and that it is stable at high temperatures only. 



CANONICAL CELL TILING MODEL 



Our Ca-Cd decoration models utilize tilings of canonical cells [18], whose vertices form 
maximal-density packings of icosahedral clusters linked along 2-fold and 3-fold icosahedral 
directions. Space is divided into 4 kinds of "canonical cells", denoted A,B,C and D. Our 
decoration rule for A,B and C cells was inferred from the refined structure of the 
Cai3Cd76.cP792 approximant [19], that we identified as the "2/1 ABC" tiling. The 
resulting A-cell decoration, when applied to the "1/1 A" tiling, reproduces the known 
MCde structure. The known quasicrystal approximants contain no D-cells and thus do not 
imply a specific decoration rule for D-cell interiors. Instead we constructed some plausible 
variants by hand and selected the lowest energy one. Our decoration rule does not impose 
specific orientations for the cluster center tetrahedra. 

The canonical cell decoration model implies decorations of prolate- and 
oblate-rhombahedra that are consistent with an earlier proposal [8] for the 1/1 cubic 
approximant. An advantage of the canonical cell approach is that it can be systematically 
extended to model higher approximants and also the icosahedral phase. 

Several canonical cell approximants were small enough that we could perform total 
energy calculations (see Fig. 2). In addition to the cubic 1/1 approximant, these were three 
rhombohedral structures: (i) "BC"-tiling (CaiaCdei), (ii) "D"-tiling (Ca2oCdio2) and (iii) 
"OR"-tihng (Ca43Cd245), where CaxCdy are formulae per primitive cell. The OR tihng 
contains all four kinds of cells (cell content AgBsCaD), and its primitive cell is a golden 
oblate rhombohedron with edge length a ~ r^Og, with "quasilattice constant" (Penrose 
rhombahedron edge length) ~5. 7A. 

ORIENTATIONAL CORRELATIONS 

We turn now to the energetics of orientational ordering in the cI176 structure. 

Consider two realizations of the cI176 structure, identical except for the orientations of the 
Cd4 tetrahedra at cluster centers. Since the tetrahedron takes two orientations, we can 
assign an Ising-like spin variable it to each one. Assign a tetrahedron the + sign if one of 
its vertices falls along the direction (1, 1, 1), and assign it a — sign if instead one of its 
vertices falls along the direction (— 1,— 1,— 1). If the tetrahedron at the unit cell vertex 
takes the + orientation and the cell center tetrahedron is also in the -|- orientation, then 
the cell vertex and center are equivalent and the symmetry is body centered, hence the 
Pearson symbol cI176. If, on the other hand, the tetrahedron at the body center takes the 
— orientation then the centering translational symmetry is broken and the stucture 
becomes primitive cubic (Pearson symbol cP176) instead of body-centered. 

Our calculations show the antiferromagnetic H — configuration is favored. To guage the 
validity of this result it is important to check convergence in the density of the fc-point 
mesh, the cutoff energy and the sensitivity to choice of pseudopotential and the density 
functional. Table I presents our study. First we vary the A;-point mesh, from 1x1x1 (the 
r point) up to 4 X 4 X 4 (all Monkhorst-Pack meshes). The table presents the convergence 
of each structure energy separately as well as the energy difference. All other 
computational parameters were held fixed during these calculations: medium precision 
(specifies cutoff energy 168 eV); ultrasoft pseudopotential; Ceperly- Alder LDA; no atomic 
relaxation. In the next series we hold constant the k-point density (we use only the F 
/c-point for speed) and test the convergence in precision going from low (cutoff 126 eV) to 
medium (cutoff 168 eV) to high (cutoff 210 eV). Next, continuing with medium precision 
and the F point, we compare unrelaxcd energies with partial relaxation (only relaxing 
clusters (a) and (b) as defined in Fig. 1) and full relaxation in which all atoms can move. 



Table I: Energy convergence studies for CaCd6.cI176 ++ and H — configurations, and their 
difference A = — All units are eV per simple cubic cell. 



Setting 


E++ 


E+- 


A 


Setting 


E++ 


E+- 


A 


1x1x1 
3x3x3 


-309.161 
-309.664 


-309.260 
-309.754 


0.099 
0.090 


2x2x2 
4x4x4 


-309.814 
-309.653 


-309.975 
-309.780 


0.161 
0.127 


low 

medium 
high 


-278.848 
-309.161 
-309.662 


-278.956 
-309.260 
-309.761 


0.108 
0.099 
0.099 


unrelaxed 

partial 

full 


-309.161 
-311.902 
-312.243 


-309.260 
-311.995 
-312.317 


0.099 
0.093 
0.074 



Table II: Energies of Ising decorations in double-length structures. Units are eV per doubled 
cell. 



Config 


Eo 


dE 


H 


Config 


Eo 


dE 


H 


+ + + + 
+ - + - 


-617.958 
-618.178 




-0.219 


C + 16 Ji + 12 J2 
C - I6J1 + 12 J2 


+ + + - 
+ + 


-618.031 
-617.986 


-0.072 
-0.028 


C + 8J2 
C + 4J2 



The maximum displacement is 0.23 A for partial relaxation and 0.30 A for full relaxation, 
always concentrated in the innermost dodecahedral shell. Specifically, those Cd atoms in 
this shell that adjoin a cluster center vacancy see the largest relaxations. 

Given the Ising-likc Z2 symmetry of the order parameter (one of two orientations) it is 
appropriate to model the energy using an Ising-model Hamiltonian. Including sufficiently 
far-neighbor interactions we can surely capture the energetics accurately. However, we have 
only a single energy difference to work with here, so wc can extract only one coupling. 
Assume this is the nearest-neighbor coupling, along the cube body diagonal, and call it Ji. 

Each BCC lattice site has 8 nearest neighbors, each of which reverses sign when going 
from -I— I- to -I — . There are two lattice sites per simple cubic cell, but we must avoid 
overcounting the bonds, since each bond is shared by two lattice sites. Hence we conclude 
that A = 16 Ji, or Ji -|-0.004 eV (using the fully relaxed value). Because the value of Ji 
is positive the interaction is antiferromagnetic. 

To determine the transition temperature for this system we wrote a simple Monte 
Carlo program to simulate the BCC Ising antiferromagnet. Actually, by simply reversing 
the sign convention for spins at body-center sites, the BCC antiferromagnet can be seen to 
be equivalent to the BCC ferromagnet. According to our simulations, the transition 
temperature should be around T=350 K. 

Because this temperature is well above the reported transition temperature of 
T=100K, we investigated the role of further neighbor interactions. If further-neighbor 
couplings have appropriate signs the magnetism can become frustrated, lowering the 
transition temperature and leading to spatial modulation of the low temperature structure. 

The next nearest neighbor lies along the cube edge, and we will call this coupling J2. 
In order to extract values of J2 we need to study a larger cell, so we doubled the cell along 
the X-axis, and considered the configurations denoted (T[(Tfo"20"| in which the cube vertex of 
the first cell has spin crj', etc. Owing to the large number of atoms present we report here 
only the results of V point calculations. 

Fitting these energies to a two-coupling Ising model yields values of Ji = -1-0.0068 and 
J2 = —0.0102 eV. Because J2 is negative, the next-nearest interaction proves ferromagnetic, 
which will increase the transition temperature and also will not lead to superlattice 
ordering. In other words, the Ising model presented so far is rather inconsistent with the 
experimental findings. 



DISCUSSION 



Given that our results disagree with experiment in both superlattice ordering and 
transition temperature it is clear that additional study is needed, especially the study of 
additional tetrahedron orientations. We restricted our attention to tetrahedron orientations 
belonging to subsets of a basic cube (based on the cI176 structure), in order to have a 
simple Ising-model description. There are many more possible orientations. Among these 
are the cuboctahedron of the cI184 structure, which our calculations (see Fig. 2) already 
show is prefered over cI176. A recent study of the MCde structure family by Gomez and 
Lidin [11] finds additional partially occupied sites, proposing the new Pearson type 
CaCd6.cI232. If an order-disorder transition exists in this new structure, it will be at a 
lower temperature than found in our initial study because: (1) the energy differences 
among orientations are generally lower; (2) the order parameter has a higher symmetry 
than the Ising spins so there is a higher orientational entropy. 
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